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Possible bound and resonant states of the hypernuclear systems Ann and AAn are sought as 
zeros of the corresponding three-body Jost functions calculated within the framework of the 
hyperspherical approach with local two-body S-wave potentials describing the nn, An, and 
AA interactions. Very wide near-threshold resonances are found for both three-body systems. 
The positions of these resonances turned out to be sensitive to the choice of the An-potential. 
Bound Ann and AAn states only appear if the two-body potentials are multiplied by a factor of ~ 1.5. 
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I. INTRODUCTION 

The A-hyperon belongs to a wide class of particles that 
are not in abundance in this world and therefore are not 
freely available for scattering experiments. The proper- 
ties of their interaction with other particles are studied 
indirectly. For example, the most important and estab- 
lished way of studying the AN interaction consists in 
measuring and calculating the spectral properties of the 
so called A-hypernuclei (see, for example, Refs. [HQ] an d 
references therein) , which are bound states of A-particles 
inside atomic nuclei. The most convenient for this pur- 
pose are very light nuclei with A < 10. Firstly, because 
such simple systems have simple spectra with only few 
well separated levels, and secondly, because they allow a 
reliable theoretical modelling based on rigorous few-body 
methods. 

The hyperon-nucleon attraction is insufficient to bind 
a AN pair. The simplest hypernucleus is therefore the 
hypertriton ^H, i.e. a bound Apn complex. Its binding 
energy is very small (the A particle is separated at ~ 
0.15 MeV) d 0, HI- So it looks like a deuteron core 
surrounded by a A- halo [1, @ . 

Similarly to traditional (non-strange) nuclear physics, 
where the deuteron is the first testing ground for any NN 
potential, the system ANN is used to constrain new mod- 
els of the hyperon-nucleon interaction. This system was 
recently analyzed in Refs. [f| @ using rigorous three- 
body equations with the potentials constructed within 
the constituent quark model. The authors of Ref. [|| 
gave another proof that the coupling between the ANN 
and TtNN channels is very important for the hypertri- 
ton binding and showed that is the only bound state 
of the ANN system. Their comprehensive analysis lacks 
only one thing: they did not consider possible three-body 
resonances. Meanwhile their results give a strong indi- 



cation that such resonances may exist and be located 
not far from the threshold energy. Indeed, they found 
that the channel Ann is attractive but not sufficient to 
produce a bound state, and the curve for its Fredholm 
determinant turnes towards zero near the threshold en- 
ergy (see Fig. 4 of Ref. @). In our present paper, we 
partly fill in the gap by considering the Ann resonance 
state. 

The AN and AA potentials are usually constructed in 
such a way that the calculations with these potentials 
reproduce experimentally known bound states of the hy- 
pernuclei. Unfortunately, it is very difficult to do scat- 
tering experiments with the A-particles because of their 
short lifetime (~ 10~ 10 sec) and extremely low intensity 
of the beams that can be obtained. 

It is well known that even when scattering data are 
available in full, it is impossible to construct an interac- 
tion potential in a unique way. One can always obtain 
different but phase-equivalent potentials (see, for exam- 
ple, Ref. In this respect the AA-case is beyond any 
hope since only few experimental points for the Ap scat- 
tering are available[8, 9]. During the decades of studying 
the hypernuclei many features of the AA-interaction have 
been revealed. However the comparison of the theoretical 
and experimental spectra remains inconclusive. Different 
potentials lead to almost the same spectra of the hyper- 
nuclei. We therefore need an additional tool for testing 
the potentials. 

In principle, such a tool could be based on studying 
the A- nucleus resonances, if they do exist [13, El- In- 
deed, while the scattering and bound states mostly reflect 
the on-shell properties of the interaction, the resonances 
strongly depend on its off-shell characteristics, which may 
be different for phase-equivalent potentials. 

Our present work is an attempt to attract the atten- 
tion of both theoreticians and experimentalists to the 
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low-energy resonances in the A-nuclear systems. As an 
example, we consider the three-body systems Ann and 
AAn in the minimal approximation, [L] = [L m i n ], of the 
hyperspherical harmonics approach. By locating the S- 
matrix poles on the second (unphysical) sheet of the com- 
plex energy surface, we show that these systems have 
near-threshold resonant states. The position of the poles 
turnes out to be strongly dependent on the choice of the 
AiV-potential. This fact supports the idea that the study- 
ing of the A-nucleus resonances could be very important 
for finding an adequate AiV-potential. 

The demands for adequate hyperon-nucleon (YN) and 
hyperon-hyperon (YY) potentials come not only from 
nuclear physics itself, but also from astrophysics. The 
studies of the neutron stars show that these very dense 
and compact objects are in fact "giant hypernuclei" (see, 
for example, Ref. [l2| and references therein). The A- 
particles appear inside neutron stars when the density 
becomes approximately two times higher than the ordi- 
nary nuclear density. The equation of state, describing 
a neutron star, involves all the inter-particle potentials 
and therefore its solutions depend on their properties. 
In particular, the strength of the short-range repulsion 
in the pairs AN and AA is crucial for determining the 
maximum mass and size of a neutron star. The repul- 
sive nature of the Ann three-body force (if it is indeed 
repulsive) would lead to additional stability of neutron 
stars. Moreover, the two-body YY interactions regulate 
the cooling behaviour of massive neutron stars [12j | . 

So, the studies of hypernuclear systems are not only 
important for reaching a better understanding of the 
physics of strange particles, but may also have an im- 
portant impact on some other branches of science. This 
is why the research in this field is carried on by many 
theoretical groups and experimental laboratories. 



only possible collision process for them is therefore the 
3^3 scattering. The wave functions describing the 
systems that cannot form clusters behave asymptotically 
as linear combinations of the incoming and outgoing hy- 
perspherical waves (see, for example, Ref. |l5[)- Thus 
it is convenient to describe the three-body configuration 
using the hyperspherical coordinates, among which only 
one (the hyperradius) runs from zero to infinity while all 
the others (the hyperangles) vary within finite ranges. 

Within the hyperspherical approach, the wave function 
is expanded in an infinite series over the hyperspherical 
harmonics (similarly to the partial wave decomposition 
in the two-body problem), and we end up with an in- 
finite system of coupled hyperradial equations, which is 
truncated in practical calculations. All the details of the 
hyperspherical approach can be found, for example, in 
the review by M. Fabre de la Ripelle [la ]. 

It should be noted that although the two-body poten- 
tials and masses for the three-body systems Ann and 
AAn are different, they can be treated using exactly the 
same equations. Indeed, in both of these systems, we 
have two identical particles with spin 1/2 and a third 
particle of the same spin. In what follows, we therefore 
consider a general system of this type. 

Let m% be the mass of one of the identical particles, 
and m,2 be the mass of the third particle. Then the total 
mass of the system is M = 2m\ + 1112 and the reduced 
masses for the identical pair and for the third particle 
are \x\ = m\/2 and ui = 2m,\mijM , respectively. With 
the Jacobi coordinates shown in Fig. [TJ the three-body 
Schrodinger equation can be written as 

( d 2 + -d r - \c? + k 2 - v) *W - (f l5 f 2 ) = , (1) 

where 



II. THREE-BODY JOST FUNCTION 



V = 2M(U n + U 13 + U 23 ) 



(2) 



There are several different ways of locating quantum 
resonances. The most adequate are the methods based 
on the rigorous definition of resonances as the S-matrix 
poles at complex energies. This definition is universal 
and applicable to the systems involving more than just 
two colliding particles. Of course, the problem of locating 
the iS-matrix poles is not an easy task, and especially for 
few-body systems. There are different approaches to this 
problem. To the best of our knowledge, so far only one of 
them has been applied to study the hyperon- nucleus reso- 
nant states. This was done in Ref. [lu] using an analytic 
continuation of the rigorous three-body equations pro- 
posed by Alt, Grassberger, and Sandhas [13| and known 
as the AGS-equations. In our present paper, we follow a 
different approach based on direct calculation of the Jost 
function using the method suggested in Ref. [14| . 

The three-body systems we consider in the present pa- 
per, namely, Ann and AAn, do not have bound states 
in any of the two-body subsystems nn, An, or AA. The 



is the sum of the two-body potentials E/y, the vectors 

{k\,k2\ represent the incident momenta of the three- 
body collision along the corresponding configuration vec- 
tors {rijfa}, the superscript [s] = ((siS2)si2S3)sa de- 
notes the spin quantum numbers for the spin-addition 
scheme s = (si + s$) + S3, the variable 
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is the hyperradius that gives the "collective" size of the 
system, k is related to the total energy, k 2 = 2ME, and 
can be called the hypermomentum, and the operator C 2 
absorbs all the angular variables. It is defined as 



f-' = -J^-4cot(2a)#- 



da 



' -4 + 



1 



sin a 



1 4 ( 4 ) 



with a = arctan(r2/ri), < a < tt/2, and £p i being 
the operators of the angular momenta associated with 
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the corresponding Jacobi coordinates. The solutions of 
the eigenvalue problem 



£ 2 Y [L] (u>) = L(L + 4)Y [l] (u>) 



(5) 



cobi vectors and the total angular momentum £ together 
with its third component m. Combining Yj^j (us) with 
the spin states X[s] = \((siS2)si2S3)s<j), we obtain the 
functions 



are the so called hyperspherical harmonics that depend 
on the hyperangles us — , flp 2 , a} including the spher- 
ical angles fl^ of the vectors f% and the angle a that de- 
termines the ratio r^/n- The subscript [L] is the multi- 
index [L] = {L,t\,t2,i,m} that includes the grand or- 
bital quantum number, 



L = t x 



2n 



0,1,2, 



(6) 



as well as the angular momenta associated with the Ja- 



( w ) = 22{ ir naa\jj a )Y lL] (cj)x[ a] 



(7) 



that constitute a full ortho-normal set of states with a 
given total angular momentum j in the spin-angular sub- 
space. 

Similarly to the two-body partial wave decomposition, 
we can expand a solution of Eq. ([1]) in the infinite series 
over the hyperspherical harmonics, 
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fei,fe 2 



r 5/2 



(8) 
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where the hyperangle sets usp and us^ are associated with 

the pairs {rj, ri} and {fci, fe}, respectively. After substi- 
tuting this expansion into Eq. (TT]) and doing the projec- 
tion onto the functions , we end up with the following 
system of hyperradial equations 



,2 , ,2 A(A+1) 



U [L][L'] = ^ V [L)[L"]U[L"][L'] , 
[£"] 

(9) 

where for the sake of simplicity we dropped the super- 
scripts jj z (indicating the conserving total angular mo- 
mentum). In Eq. @, 



and A = L + 3/2. Since we consider a system that can- 
not form clusters, the asymptotic behaviour of its wave 
function may only involve the incoming and outgoing hy- 
perspherical waves ~ exp(=pz/cr), which are the products 
of the corresponding spherical waves along the Jacobi 
radii r\ and r^. 



e *iri e ik 2 r 2 



gikr cos ot^ikr sin a ^ zr 



V {L]W] (r)=2M 
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t- U 23 ) (w)da; 

(10) 



We therefore look for the solution of matrix equation © 
as 



u [L] [L/] (E, r) = (kr)F^> LI] (E, r) + h?> (*r)^] r) 



.(in) 



(+), 



p(out) 



(11) 



where the incoming and outgoing hyperspherical waves 
described by the Riccati-Hankel functions, 



hf\kr) 



\kr\ 



=Fz exp [±i{kr — Xtt/2)] , 



(12) 



are included explicitly. The matrices Fmiiffl (E, r) are 
new unknown functions. In the theory of ordinary differ- 
ential equations, this way of finding solution is known as 
the variation parameters method (see, for example, Ref. 



Since instead of one unknown matrix 



troduce two unknown matrices F, 



(in/out) 



l[L][L>] 



we m- 



'[£][£'] ' ^ ne y cann °t be 
independent. We therefore can impose an arbitrary con- 
dition that relates them to each other. As such condition, 
it is convenient to choose the following equation 



(kr)d r F$ [L1] (E, r) + h™ (kr)d r F^ (E, r) 



(+) 



r (out) 



= 



(13) 
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which is standard in the variation parameters method the condition (fT"3] ). we obtain the following system of first 
and is called the Lagrange condition. Substituting the order equations for these unknown matrices 
ansatz (fTTj) into the hyperradial equation ([9]) and using 



°r r [L][L' 



a p(out) _ 
°r*[L][L'] - 



,(+) 



V, 



[L][L"] 



2ik ^ 

[L»] 

h (-) 

Z2 V [L][L"] 
[£"] 



2ik 



,l X» e \L"\\L> 



n \" *[L»][L>] 



, (+) p (out) 
a X" r [L"][L> 



(14) 



which are equivalent to the second order Eq. @. The 
regularity of a physical wave function at r = implies 
the following boundary conditions 



F, 



(in) 
[L][L 



, ] (E 7 0) = F^, ] (E,0) = S [L][L/] . (15) 



With 



these conditions, the columns of the matrix 
U [L][L'] (E, r) are not only regular but linearly indepen- 
dent as well. Therefore any regular column <j>nn (E, r) 
obeying Eq. ([9]), can be written as a linear combination 
of the columns of matrix U[M[£/i (E, r). In other words, 
the matrix u^^iAE^r) is a complete basis for the regu- 
lar solutions. 

At large hyperradius where the potentials vanish, i.e. 

V [LW ](r) — 0, (16) 



the right-hand sides of Eqs. ([M]) should tend to zero 
and therefore the matrices F^^^\e, r) converge to the 
energy-dependent constants, 



f (in/out)/ m 
J[L][L>] KP) 



lun F$(°^(E,r), 



(17) 



that by analogy with the two-body case can be called the 
Jost matrices. The convergency of these limits, however, 
depends on the choice of the energy E and on how fast 
the potential matrix Vrnrx,/](r) vanishes when r — > oo. 

When the energy is real and positive (scattering 
states), the vanishing of the right hand sides of Eqs. (fl4|) 
at large distances is completely determined by the be- 
haviour of Vrt][£/] (r). It can be shown that in such a case 



the limits (|17|) exist if V[l][l'] ( r ) vanishes faster than 1/r. 

With negative and complex energies there is a technical 
complication. The problem is that one of the Riccati- 
Hankel functions on the right hand side of Eqs. ([M]) 
is always exponentially diverging. Therefore, if at large 
distances the potential matrix vanishes not fast enough, 
the convergency of (jTTJ) is not achieved. This problem 
can be easily circumvented by using different path to the 
far-away point (see Fig. \2§ ■ 



This is known as the the complex rotation of the co- 
ordinate. All the details concerning convergency of the 
limits (|17|) and the use of c omp lex rotation for this pur- 
pose can be found in Refs. H El |H El HI • 

As was said before, the columns of the matrix function 
U [L][L']{E , r) constitute the regular basis using which we 
can construct a physical solution </>rn (E, r) with given 
boundary conditions at infinity, 



4>[L](E,r) = ^2u [L][L >](E,r)C [L r 
[£'] 



(18) 



where C[£] are the combination coefficients. 

The spectral points E n (bound and resonant states) are 
those at which the physical solution has only outgoing 
waves in its asymptotics, i.e. when 



(19) 



This homogeneous system has a non-trivial solution if 
and only if 



det/K L ,](£„) = 0, 



(20) 



which determines the spectral energies E n . As can be 
easily shown 23], the S'-matrix is given by 



f( in )/ 



and therefore at the energies E n it has poles. 



(21) 



III. TWO-BODY POTENTIALS 

In our calculations, we used local two-body potentials 
describing the interaction between two neutrons, A and 
neutron, and between two A-particles. For all these po- 
tentials, we used the same functional form, namely, 
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U(p) 



1 + P a 1 - P a 

Mp) t^Mp) t^Mp) 



A n (p) = W n exp(-a n p 2 ) , n = l,2,3, 



(22) 
(23) 



r 



where P a and P r are the permutation operators in the 
spin and configuration spaces, respectively. The form 
of U(p) as well as the parameters were taken from Ref. 
[24[. In order to explore how sensitive the positions of 
the three-body resonances are to the choice of underlying 
two-body potentials, we did the calculations with three 
different sets of parameters for the An-potential. All the 
sets of parameters we used, are given in Table HI 



IV. THE MINIMAL APPROXIMATION 



The system (|T4|) consists of infinite number of equa- 
tions. For any practical calculation, one has to truncate 
it somewhere. Before going any further, it is very logical 
to try the simplest approximation, namely, when only 
the first terms of the sums on the right hand sides of 
Eqs. are retained. This corresponds to the mini- 

mal (n = 0) value of the grand orbital number © and is 
called the hypercentral approximation, [L] = [L m in]- We 
assume that the two-body subsystems are in the S'-wave 
states (£i =£2 = 0), which means that 

A = A m i n = — . 

So, in the minimal approximation, instead of the infinite 
system (j 14[) . we remain with only one equation, 

i{E,r) = 2M{U)u(E,r) , 
(24) 

where all unnecessary subscripts are dropped, and the 
brackets on the right hand side mean the following inte- 
gration 



q2 _|_ ^2 A m ; n (A m i n + 1) 



*[L*JL]M (^12 + ^13 + ^23) fcf^Hdw • 

(25) 



(U)(r) 



From the mathematical point of view, Eq. (|24|) looks 
exactly like the two-body radial Schrodinger equation. 
The only difference is that the angular momentum is not 
an integer number. 

The explicit expression for the integral (|25[) is given in 
the Appendix. The hypercentral potentials (U) for the 
systems Ann and AAn are shown in Figs [3] and 2] With 
these hyperradial potentials the corresponding differen- 
tial equations determining the three-body Jost functions, 
were numerically solved with complex values of the en- 
ergy. The results of these calculations are discussed next. 



V. NUMERICAL RESULTS AND CONCLUSION 

When looking for zeros of the three-body Jost func- 
tions, we found that there were no such zeros at real neg- 
ative energies. In other words, neither the system Ann 
nor AAn have bound states. 

The only zeros we found were located on the unphysi- 
cal sheet of the energy surface, in the resonance domain. 
The resonance energies are given in Tables [TTI and Mil and 
shown in Fig. [5l As is seen, the positions of these reso- 
nances depend on the choice of the An potential. For the 
choice "C" , the resonances become sub-threshold. 

In order to estimate how far our three-body systems 
are from being bound, we artificially increased the depths 
of the potentials by multiplying them by a scaling factor. 
When this factor was increased from 1 upwards, the Jost 
function zeros moved towards the origin of the energy 
surface. At the value of approximately 1.5, the zeros 
crossed the threshold and moved onto the real negative 
axis. In other words, the bound states can appear if the 
potential strength is increased by ~ 50%. 

The fact that we did not find bound Ann or AAn states 
is not surprising at all. As is shown in Refs. [!, Q , the 
system ANN in the state with the three-body isospin 1 
and spin s = 1/2 is not bound even when the virtual 
processes of A — S conversion are taken into account, 
although this conversion increases the attraction in the 
system. Simple but convincing argumentation of Ref. 
[251 ] leads us to the conclusion that the AAn system also 
cannot be bound. Indeed, the system AAn is a "mirror" 
image of Ann, where the A and n replace each other. This 
means that the potential term U = U nn + U\ n + Uau of 
the three-body Hamiltonian is replaced with U = Uaa + 
U\ n + U\ n . Since the attraction of Uaa is weaker than 
that of U nn , we may conclude that the system AAn has 
less chances to be bound than the system Ann. The 
calculations performed in Refs. 0, Hfl H3, [Hj], show 
that even the heavier hypernucleus a/J^ (i-e. the system 
AApn) is bound very weakly, if bound at all. 

Multiplying the two-body potentials by an appropri- 
ate scaling factor, we can always generate an artificial 
three-body bound state, i.e. a pole of the S'-matrix on 
the physical sheet of the i?-surface at a negative energy. 
Apparently this pole cannot disappear when the scaling 
factor returns to its natural value of 1. The pole sim- 
ply moves via the threshold onto the unphysical sheet. 
Since both the systems we consider, are not far from be- 
ing bound, their corresponding poles cannot be far away 
from the threshold energy. And indeed we located them 
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at low energies. 

What we found is, of course, an estimate. But it clearly 
shows that there are near-threshold resonances of the sys- 
tems Ann and AAn. Actual location of the poles most 
probably is more close to the threshold energy. An in- 
clusion of the channels AA~ — UN and AA — SAT would 
definitely increase the attraction in our systems (see Ref. 

) and this would make the widths of the resonances 
smaller. 

As we have demonstrated, the positions of the reso- 
nances strongly depend on the choice of the two-body po- 
tentials. If such resonances are observed experimentally, 
they may serve as an additional instrument for construct- 
ing adequate YN and YY potentials. There are many 
possible reactions where the three-body resonances Ann 
and AAn may manifest themselves. As an example, we 
can mention the inelastic collision of the K~ meson with 
the a particle, 

K~ + 4 He — > p + A + n + n, (26) 

that produces a proton and the system we are looking for. 



If a short-lived cluster Ann is formed in the final state 
of this collision, it should be seen in the corresponding 
two-body kinematics p — Ann. The processes of the type 
(|2"6")) fall under the experimental programme AMADEUS 
[30| (in the INFN, Italy) and, in principle, this reaction 
could be thoroughly studied. 



APPENDIX 

Hypercentral potential (|25|) consists of the three terms 
(U) = (t/12) + (C/13) + (U23) , (27) 

where Uij is the two-body potential acting between par- 
ticles i and j. As was mentioned above, we can consider 
both nnA and AAn systems in a unified way. Let 1 and 

2 be the identical particles, i.e. the nn or AA pair, and 

3 be the remaining A-particle or neutron, respectively. 
The six-dimensional volume element is 



dr\dr2 = r 2 r 2 dr idr2 sin 9 f 1 d9p 1 d(pp 1 sin 9 f 2 d9p 2 d<ff 2 

= r 5 dr— sin 2 (2a)da sin 9? t d9? 1 d(p? 1 sin 9p 2 d9? 2 dip?, 
= r 5 drdoj . 

I 



Therefore in the five-dimensional integral (|25| the volume 
element is 

du> = — sin 2 (2a) sin 9^ sin 9? 2 dad9? 1 d<pp 1 d9j? 2 dtpF 2 . 

' (28) 
Since we assume that l\ = £2 = and L = L mnl = 0, 
the sum (J7J is reduced to a single term, 

$g; in] ( W )=y [imin] (a;)x [s] , (29) 

where the quantum numbers jj z coincide with sa. The 
two-body spin S12 of the identical pair in the S'-wave state 
must be zero. As a result the three-body spin s is always 
1/2. The hyperspherical harmonics Yji m . n ]([j) is trivial 
(independent of the angles), 

Y- [Lmin] H=7r- 3 / 2 , (30) 



which means that the action of the permutation opera- 
tors P r for all three terms in Eq. (|2T|) is also trivial: its 
eigenvalue is 1, 



P[jY [Lmin] = Y [Lmia] , ij = {12}, {13}, {23} . (31) 
The spin permutation operator P° 2 for the identical pair 
{12} changes the sign of X[s\i 



P?2X[s] - -X[s] , (32) 



because s\2 = in [s] — ((siS2)si2S3)s<7. For the other 
two pairs, its action is a bit more complicated. Indeed, 
recoupling the spins, 



\((sis 2 )si2S3)scr) = ^2 \((s3Si)s 3 is 2 )scr)(((s3Si)s3iS2)sa\((siS2)si2S 3 )scr) 
S31 
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= |((s 3 Si)Os 2 )scr) 



1 i 
? ? J 

2 2 u 
1 



v / 3|((s3Si)ls 2 )so') 



2 2 



|((s 3 Si)0s 2 )scr) 



|((s 3 Si)ls 2 )so-) 



When inserting the potentials E/jy given by Eq. (|22[) . into 
the integral (|23|) , we should use the following interparticle 
(33) distances (see Fig. [1]), 

Similarly, it is easy to find for the pair {23} that 

xt^Xls] = \ ■ (34) 



we find that 



Pl2 




(35) 



P13 



M . 2 M 

— sin a H cos- 2 a 

M2 



M 



2VMiM2 



sin(2a) cos 0f? 2 



(36) 



P2 3 



/'2 



sin 2 a 



M 
4^i 



ill 



2^/MiM2 



sin(2a) cos 0p 2 



(37) 



r 



Since the particles 1 and 2 are identical the interactions 
J7i3 and t/ 23 are the same. Moreover, according to Eqs. 
([3"3"Il3~i|) , the product xf s ]U~i3X[s] has the same dependence 
on pi 3 as Xm^23X[ s ] depends on /9 23 . Actually, they de- 
pend on p 2 3 and p 23 because all the terms in the potential 
(|22| are of the Gaussian form. Comparing Eqs. (|36|) and 



(|3"T)) . we see that the integrands for (U13) and (C/ 23 ) differ 
only in the sign of the power of the exponential factors 
corresponding to the last terms of p 2 3 and p\ 3 . This dif- 
ference however has no effect on the integrals. Indeed, 
the integration over 8? 2 , 



J exp(±/cos#f 2 ) sm9? 2 d9? 2 — J exp(±ft)dt 



-f\ - 



= J sinh(/) , 



gives the same result for both signs. Therefore (U13} — (C/ 23 ) and hence 

(U) = (C/i2> + 2(U 13 ) . 



(38) 



Performing trivial integrations over (p^, tpp 2 , 9p 1 , and 9? 2 (trivial in the case of (U12}), we obtain the following 
expressions for the terms of the hypercentral potential 



(U 12 ) = 1 f ' dasin 2 (2a) [v^ {12} exp (-aj 12 V 2 ) - V^ 3 {12} exp (-a^ 12} 



rjr 



(39) 
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(Uia) = ^£ dasin 2 (2a)[iy i {13} exp(-al 13} Cr 2 ) S (al 13} er 2 ) (40) 

-^f 3 >exp(-4 13 >C, 2 ) S (4 13} ^ 2 ) 
^3 {13} exp(-4 13 >Cr 2 ) S (4 13 >^^ 



4 

I 

where The parameters and , where the symbol {ij} 

j. means a choice of the pair of interacting particles, are 

^(qA = — cos 2 a ] given in Table Q] For each (complex) value of the hy- 

^ perradius r, which was needed in our calculations, the 

integrals p9|) and (|40)l were evaluated numerically. 

>/ \ • 2 2 

uoO = — sin a + - — cos a, 

M2 4/il 
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Tin 


AA 


An (A) 


An (B) 


An(C) 


Wi (MeV) 


200.0 


200.0 


200.0 


600.0 


5000 


W 2 (MeV) 


178.0 





106.5 


52.61 


47.87 


Wi (MeV) 


91.85 


130.8 


118.65 


66.22 


61.66 


ai (frn^ 2 ) 


1.487 


2.776 


1.638 


5.824 


18.04 


a 2 (fm- 2 ) 


0.639 





0.7864 


0.6582 


0.6399 


a 3 (fm" 2 ) 


0.465 


1.062 


0.7513 


0.6460 


0.6325 





1 


1 


1 


1 


1 



TABLE I: Parameters of the potential (|22[1 for the pairs nn, AA, and An. For the system An, three different sets of parameters 
(denoted as A, B, and C) are given. All the parameters are taken from Ref. [24l ]. 



An-potential 


A 


B 


C 


E , (MeV) 


0.551 - ^4.698 


0.456 - ^4.885 


-0.149 - ^5.783 



TABLE II: Complex resonance energies Eq — E T T for the system Ann with the three choices of An-potential. 
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An-potential 


A 


B 


C 


E , (MeV) 


0.096 - ^8.392 


0.034 - ^8.438 


-0.552 - ^8.681 



TABLE III: Complex resonance energies Eo = E T — —T for the system AAn with the three choices of An-potential. 
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FIG. 1: Jacobi vectors defining the spatial configuration of a three-body system of two identical (filled circles) and one different 
(open circle) particles. 



Im r 







Rl 











R 



Rer 



FIG. 2: Deformed contour for integrating Eqs. (|14p from r = to r — R when the energy is complex. 
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FIG. 3: The hypercentral potential given by Eq. I|25p for the system Ann with the three choices (A, B, and C) of the An interaction. 




FIG. 4: The hypercentral 



potential given by Eq. I|25p for the system AA71 with the three choices (A, B, and C) of the An interaction. 
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FIG. 5: Resonance points for the systems Ann and AAn found on the unphysical sheet of the energy surface with the three sets (A, 
B, and C) of parameters of the An-potential given in Table© 



